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We suggest certain effects, caused by interaction between rotation and gravita¬ 
tion with turbulence structure, for the cooling/heating of dispersed phase of non- 
isothermal particles in rotating turbulent fluid flows. These effects are obtained 
through the derivation of kinetic or probability density function based macroscopic 
equations for the particles. In doing so, for one-way temperature coupling, we also 
show that homogeneous, isotropic non-isothermal fluid turbulence does not influence 
the mean temperature (though it influences mean velocity) of the dispersed phase of 
particles settling due to gravitational force in the isotropic turbulence. 


I. INTRODUCTION 


Different phenomena of non-isothermal two-phase flows existing in situations, such as, 
dispersion of aerosols, dust devils, dust storms and their effects on atmosphere, protoplane¬ 
tary nebula and disks, birth and growth of cloud droplets, all involve rotation, gravitational 
held (e.g. see [l-5|). Researchers continue to face challenges in accurately explaining the 


observed and discovering the new phenomena 


111 ] due to the involved turbulence closure 


problem with complexity further enhanced by the presence of dispersed phase of particles 


and droplets. In recent years, com 
huid turbulence closure problem 


Dination of kinetic approach 


12j and theories for classical 


13| has evolved in a mathematical rigorous approach of 


probability density function (PDF) capable of capturing certain phenomena for the dispersed 
phase of particles/droplets in turbulent hows (e.g. see reference llj] and references cited 
therein). In our ongoing ehorts toward unihcation of various aspects and explanation of the 
phenomena of two-phase turbulence [w, 11, 14|, in this letter we apply PDF approach 
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2l| to important situation of two-phase rotating, non-isothermal, turbulent flows yielding 


new effects for collective behavior of particles’ temperature. 


II. ANALYSIS 


Consider particles dispersed in non-isothermal rotating flows moving due to forces of 
fluid drag, rotation (Coriolis and centrifugal) and gravity. Particles also exchange heat with 
the surrounding fluid and are considered here as point particles for analysis purpose. The 
trajectory and temperature of each particle is governed by the Lagrangian equations for its 
position X, velocity V and temperature Tp, written in a coordinate system having constant 
angular velocity f2 as 


dXi 

dt 




d^ 

dt 


— Pv{Ui — Vi) + 2eiabVaV,i, + fi, 


dTp 

dt 


= fie{T-Tp) + Q{Tp). 


( 1 ) 


( 2 ) 


Here, = l/r^, Tp is particle velocity time constant, Ui and T are carrier fluid velocity and 
temperature in the vicinity of the particle, eiab is the Levi-Civita’s alternating tensor, and 
fi = ff -|- ff account for components of centrifugal V = —Li x (il x r) and gravitational 
f 3 _ jgj. accelerations with = —{GM/r^) and position vector of particle r and r =| r |. 
Also, G is the gravitational constant, M is the mass responsible for f^, is inverse of 
the particle temperature time constant, Q{Tp) is a function of Tp accounting for particle 
temperature source and/or heating and cooling of particles due to radiation, e.g. see cases 
of optically thick and thin limits in situation protoplanetary disk jd]. The collective behavior 
of particles can be represented by the Eulerian instantaneous equations in physical space x 
and time t for particles number density ?t,(x, t), velocity f/(x, t) (obtained from Eq. ([1]) Q) 
and for temperature 0 which can be formed from Eq. ([2]) in an usual manner, written as 


dQ de 

— + v,— = fdg{T-e) + Q{e). 


dt 


dxi 


( 3 ) 


To extract statistical information for collective behavior of particles from the instantaneous 
equations, the density-weighted average (denoted by overbar) of the equations is preferable 
instead of non-weighted average due to the established reason ll|. The mean number 
density N = {n) then satishes ON/dt + d{NVj)/dxj = 0, density-weighted averages Vi and 
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0 satisfy 


DtV, 


Qv' v'- _ _ 

= /3,{{U,) - V,) + Rj + 


DtQ + ^ = M{T) - 0) + Q - v'J’ 


d\nN — 

-7, -+ Pet" 

ox„. 


(4) 


(5) 


with Rj = 2ejabVaP^b + fj — v'^v'jdlnN / dxn and operator Dt = d/dt + Vidjdxi. Also for 
any instantaneons variable A, A = {nA)lN with ensemble average denoted by (), a' and 
a" are flnctnations in A = A + a' = (A) + a" over A and (A), respectively. In Eqs. (jl]) 
and ([5]), u" and t" represent respectively, flnctnations in flnid velocity and temperatnre, 
in the vicinity of particle, over the flnid mean velocity (Ui) and temperatnre (T). The 

from the 


appeared nnknown terms v'^v'j^v'^O' and, in particnlar, u"R" are difficnlt to mode 


instantaneons eqnations bnt can be obtained with ease in PDF approach 
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M 


llTj. 


III. NON-ISOTHERMAL, ISOTROPIC TURBULENCE CASE 


Before obtaining expressions for the nnknown terms, consider Eqs. (jl]) and ([5]) for an 
ideal sitnation homogeneons, isotropic, non-isothermal flnid tnrbnlence with nniform (T), 
Ob = 0 and having dispersed particles settling nnder the gravity with Q = 0 , /* = — 
where g is acceleration dne to gravity. The eqnations simplify to U 3 = —g/Pv + u'l and 
0 = (T) +t" having additional drift velocity u'l and temperatnre f", for the particles, cansed 
by the tnrbnlence velocity and temperatnre structures. For the isotropic case, though u'l 
is nonzero 22 |, we show now that t" = 0 when the effects of particle temperature on the 
fluid temperature are neglected (i.e. one way coupling of temperature). For the one-way 


temperature coupling and diffusivity constant a, t" is governed by ^ + 


= a-. 


and 


which suggests that for each realization of u", at any time f, t" and —t" are equally probable 
at any location of the dispersed particle. This suggests that the density-weighted average 
of t” vanishes i.e. t" = 0. All these imply that particles’ mean temperature 0 = (T) is 
identical in both the cases of isotropic non-isothermal turbulent and stationary fluids, both 
having identical uniform mean fluid temperature (T). 
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IV. PDF APPROACH 


We use PDF approach to obtain expressions for unknown terms and, in particular, u" 

and t" in general situation of turbulent flows. The equation for ensemble average of phase 

space density W (x, v, 0, t) of the particles, representing PDF equation, is 

d(W) d d d 

at + + fd{W) 


~[/?„({r) - 9) + Q\(W) = -±[fi,(u';w)\ - ^IA(W)], 


( 6 ) 

which can be obtained from Eqs. o and (El) and Liouville’s theorem [l^. Here x, v, and 
9 are phase space variables corresponding to X, V, and Tp, respectively. The expressions 
for unknown terms {u'lW) and {t''W) appearing in Eg. dp c an be obtained by employing 


Furutsu-Donsker-Novikov functional formula 


18 


, |20|, |23|, [2^, which are exact expressions 


when u'l and t" along the particle path have Gaussian distribution. The expressions are 
^{u'lW) = - 


d d d 

v;- ^ki + T;- f'ki + 

uX]^ uVj^ Uu 


li 


Mfw) = 

where various tensors are 


d d d 

— —A^ T ——H/j. T ——D 
dxk dvk 09 


w, 

(W). 


(7) 

( 8 ) 


dt2{t"Uj{t\t2))Gjk{t2\t), 


^ dt2{t"u'-{t\t2))^Gjk{t2\t), 


fffc l^vl 

r = /3y/30 dt 2 (u''{t\t 2 )dt''/dxk) Gjk{t 2 \t) 

D = /dy/de f dt2{t''u''{t\t2))Gj{t2\t) + Q,2 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 


with D 2 = Idg dt 2 {t"t''{t\t 2 ))G^{t 2 \t) and expressions for remaining tensors can be obtained 
by writing jdyu'l instead of jdot” on the right-hand side (rhs) of Eqs. (l9|)- (fT^ and in D 2 and 
changing A^. —)■ A^i, H^ —)■ D ^ a;* and F —)■ 7 *. In these expressions, shorthand no¬ 
tations u'l{t\t 2 ),t''{t\t 2 ) represent m"(x, v, 6*, t|t 2 ), ^Gl^ 2 ), the argument (x, v, 6*, t|t 2 ) 

represents the value of u'l and t" at time ^2 in the vicinity of particle that passes through x 
at time t with velocity v and temperature 9. Also, u" and t" represent ul{'x,t) and t"(x, t), 
respectively and now onwards too. The equations for Gjk(t 2 \t), Gj(t 2 \t), G^(t 2 \t) V t > ^2 are 


'DGjk{t2\t) — (dyG 




d{Uk) 

Oxi 


dGja 

j, 

at 


Ofk 

OXi 


Gji — 5jkS{t — ^ 2 ) 


(13) 
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j^G\h\t) + PeG^ - = 5{t - t2) 


( 14 ) 


(15) 


where operator "D = cP/dt^ + I3vd/dt. For later convenience, Aj and F are written in different 
forms as 

^ = [ ds{t"Avi), ^ = [ ds{Avidt"/dxi) (16) 

Pe JO Pe JO 

where An* = Jq dt 2 represents change in particle velocity due to /5^M^'(t|t2), 

during time 0 to s, along the trajectory that passes through x, v at t. 

The various tensors contain statistical properties related to u'' and t" along the particle 
oath and which can be obtained from the modelled equation for u” along the particle path 
ll| . exact equation for t” 


dt 


d{T) 


V72j.// // 

—-aVf 


,dt 


dt" 




dxi 


(17) 


with model aVH” = —P'/Tq based on the simple lEM model [25| suggested for single-phase 
turbulence. It should be noted that integral time scale Tg is different than the time scale of 
lEM model. To obtain {t''u'-{t\t 2 ))i multiply Eq. flTTll by u”{t\t 2 ) and take ensemble average 
to yield 1 ^ 2 ))/dt = -{fu”{t\t 2 ))lfe +Sj wit\i Sj = -{u”{t\t 2 )[u'l^+ {Ui-Vi)^l\). 

Further, the first order solution for {t''u”{t\t 2 )) can be obtained by neglecting Sj and the 
result is an exponential form {t''u”{t\t 2 )) = {f {t\t 2 )u'-{t\t 2 ))e~P~^'^P^^ which is considered 
here for further analysis. 

The continuum equations for the dispersed phase can be obtained by taking various 
moments of Eq. ([6]) after substituting for P^{u[W) and from Eqs. ([7]) and ([H]). The 

equations for E, = A j Vi{W)dvd9 and 0 = A j 9(pv)dvd9 are identical to Eqs. dl]) and 
o. respectively, with 

_ T... A/V f) _ 

(18) 


\jdN d- 


Pet" = 


AkdN 


N dxk dx. 


dAi T. 
+ F. 


(19) 


Though the equations for higher order correlations v'^Vj,v!^9' etc. can also be derived from 
the PDF equation, we focus our discussion on different phenomena contained in Eq. flT^ 
for particles’ temperature. 
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V. NEW EFFECTS 


The phenomena contained in various terms of Eq. fflSD for additional drift of particles 
are described earlier [h|. Here we discuss phenomena related to temperature of collective 
particles as described by Eq. ([T9]). For the discussion, we consider density weighted average 
of various tensors to be equal to their respective instantaneous values, e.g. Xij = \ij. It 
should be noted that the form of Eqs. dSj) and flT^ remains unchanged irrespective of whether 
the rotation and gravitation exist or not. The rotation and gravitation effects are felt by 
these equations through Gij and Gj appearing in various tensors. The hrst term on the rhs of 
Eq. fll9p represents phenomenon of additional increase or decrease in mean temperature of 
particles caused by combination of statistical nature of the turbulence structure represented 
by Ak and gradient of particles mean number density in general situation of non-isothermal 
two-phase turbulent flows. The rhs of Eq. (IT^ vanishes when fluid variables along the 
particle path are correlated by delta function in time as suggested by expressions for and 
T along with Gjk{t\t) = 0 [n|. When the correlation time is finite, the last two terms of Eq. 
fll9p account for certain phenomena which we discuss now for cases of (1) slow rotation i.e. 
is small and (2) fast rotation without gravity i.e. /f = 0. 


VI. SLOW ROTATION CASE 


For slow rotation, we simplify the last two terms in Eq. (fT^ by expanding the solution 
of Eq. flT^ as Gjk = G^^ + ^]k G^^ = G^{t 2 \t)Sjk = Sjk[l — and 


VG], - ^G%^^ - 


G . = 0, 


( 20 ) 


using exponential form for {t”u''{t\t 2 )) and {t"{t\t 2 )u'l{t\t 2 )) = {f'u”). The simplification 
yields 

-^Ai-hT = -I3^l3e f - (3^(3e [ dtiGl^{ti\t)e~d-^^dT9 

OXi Jo OXi Jo 


- {t"dul/dxk)l3eA 


( 21 ) 


where A = {T^)! — e [e h/^^+i/A) _ x]}. Using the solution for Gl^ from Eq. 

), the second term on the rhs of Eq. fl2Tp can be further simplihed as 

- - 2f3,f3ee,kbnMt) - (22) 
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where A\ = J^dtile dt 2 G^{t 2 \t)G^{ti\t 2 )] and A 2 = 

/Q(iti[e /j* (it 2 G°(t 2 |^)^G°(ti|t 2 )]- The first term in (1^ represents contribn- 
tion to heating/cooling of particles dne to mean shear rate. The second term in fl2^ 
snggests that correlation between t" and the component of flnid vorticity at particle 
location in the direction of rotation produces additional heating/cooling 

for particle phase and whose origin is the Coriolis force. The last term in fl 2 ^ represents 
heating/cooling of the dispersed phase due to the interaction of centrifugal and gravitational 
forces with turbulence fluctuations u" and t”. For further discussion, we consider a finite 
value for rotation only about 2:3 axis he. = 1^2 = 0 and ffs 7 ^ 0. Then, the centrifugal 
part of fk is = Xkfi^ — 5kzXk^\ and its contribution to the last term in fl 22 ]) becomes 

- l3^l30A'^{t)nl{f {dulldxk - du'l/dx'i}), (23) 

which describes the phenomenon of heating/cooling of the dispersed phase in the presence 
of a centrifugal force. The effects of gravitational force contained in the last term in fl2^ 
can be simplihed as 

- ^^PeA\{t)[{t''dul/dxk)f^ + {t''dul/dxi)x~-] (24) 

OXk 

where the last term arises due to variation of gravitational acceleration in space. 


VII. FAST ROTATION WITHOUT GRAVITY 


For large value of t, Eqs. flTB]) suggest that most contributions to Aj and F come from 
the correlations of t” and dt”/dxi with velocity of particles, between s(< t) and t, present 
in the region near to x and t and passing through x at t. Assuming exponential form 
for these correlations with integral time scale Tg, the last two terms in Eq. (IT^ become 


approximately equal to —fde{t"dVi (x, t) / dxi) /g ds e and which, along with an approximate 
2t U e I 

relation rotation without gravity [ 26 1 ], suggests phenomenon of 

additional cooling/heating of particles caused by fast rotation. This phenomenon does not 
diverge with the growth of the angular velocity. 









VIII. EFFECT OF FLUID COMPRESSIBILITY 


Fluid compressibility gives rise to interesting effects, of turbulent thermal diffusion and 
barrodiffusion of particles, which are contained in Eq. (ITSll in| . For particles dispersed in 


a compressible ideal gas, compressibility effects on the temperature of the particles can be 
obtained by using the following relation [l^ 


{t''du'l/dxk) = {t” u'l)[d\n.{T) / dxk — d\n.{P) / dxu]. 


(25) 


Substituting it into the last term in Eq. (]2l|l suggests turbulent thermal and pressure 
effects on the heating and cooling of the dispersed phase analogous to the two phenomena of 
turbulent thermal diffusion and barodiffusion Substituting Eq. (|2^ into Eqs. fl23l) and 


suggests that the first term on the rhs of these equations further adds to these thermal 
and pressure effects in the presence of rotation and gravitation held. 


IX. CONCLUDING REMARKS 


It has been our goal to accurately describe and discover the phenomena related to im¬ 
portant situation of dispersed phase of non-isothermal particles in rotating turbulent hows 
under the inhuence of gravitational held. The mean temperature 0 of the particle phase 
is found to be ahected by the term t" in Eq. ([5]) containing new ehects caused by mutual 
interactions of turbulent structure of huid temperature and velocity with the rotation and 
gravitation held. These ehects and their accurate description have been surfaced through 
the application of mathematical rigorous PDF approach having established capability in 
capturing certain phenomena (see ll| and references cited therein). 
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